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IN THE UNITED STATES ELECTED OFFICE 
OF THE UNITED STATES PATENT AND TRADEMARK OFFICE 
UNDER THE PATENT COOPERATION TREATY-CHAPTER II 

" PRELIMINARY AMENDMENT " 

APPLICANT: Stefan SCHAFFLER et al. 

SERIAL NO.: EXAMINER: 
FILING DATE: ART UNIT: 

INTERNATIONAL APPLICATION NO.: PCT/DE99/03834 
INTERNATIONAL FILING DATE: 1 December 1999 

INVENTION: METHOD AND DEVICE FOR REDUCING A NUMBER 

OF MEASURED VALUES OF A TECHNICAL SYSTEM 

Hon. Assistant Commissioner for Patents 
Box PCT 
Washington D.C. 20231 

SIR: 

Amend the above-identified international application before entry into the 
national stage before the U.S. Patent & Trademark Office under 35 U.S.C. §371 
as follows: 

IN THE SPECIFICATION 

Please substitute the specification in the file with the enclosed substitute 
specification in compliance with 37 CFR 1.125(b). Furthermore, a separate 
marked up copy of the specification that shows all changes relative to the previous 
specification is also enclosed. 
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IN THE CLAIMS 

Please cancel all claims without prejudice and add new claims 1 1-20 as 
follows. 

WE CLAIM: 

11. A method for modeling a technical process of an engineering plant 
comprising: 

measuring an initial set of empirical values at various steps of a technical 
process using sensors while said technical process is operating based on a 
predetermined set of parameters; 

screening out a set of empirical values from the initial set of empirical values 
for reducing a size of the initial set of empirical values to obtain a screened set of 
empirical values by: 

dividing the initial set of empirical values into classes based on a predefined criteria, 
followed by assessing each empirical value in each class with respect to a predefined 
first threshold value, and if a result of said assessing step lies below said predefined 
first threshold value, then screening out said empirical value, further assessing each 
class with respect to a predefined second threshold value, if a result of said 
assessing step lies below the second predefined threshold value, then, screening out 
said class; and 

modeling said technical process using said screened set of empirical values 
to obtain a model result. 

12. A method according to claim 11, wherein the predefined criteria in said 
dividing step is based on the predetermined first set of parameters. 
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13. A method according to claim 1 1, further comprising the step of : 
determining an empirical value associated with a transient phase of the 

technical process resulting from a modification of the predetermined set of 
parameters; and 

5 screening out the empirical value associated with the transient phase. 

14. A method according to claim 1 1, further comprising the step of: 
reducing a number of empirical values in a class by selecting a representative 

empirical value for the class. 

15. A method according to claim 1 4, wherein the representative empirical value 
10 is an average of the empirical values in the class. 

16. A method according to claim 1 4, wherein the representative empirical value 
is one of a maximum value and a minimum value of the empirical values in the class. 

17. A method according to claim 14, wherein the representative empirical value 
is a median of the empirical values in the class. 

15 18. A method according to claim 1 1, further comprising the step of: 

screening out a class with fewer number of empirical values than a predefined 
number. 

19. A method according to claim 1 1 , wherein the result of said assessing step is 
a difference of the empirical value in the class with the predefined first threshold 
2 0 value. 
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20. A system for modeling a technical process of an engineering plant 
comprising: 

a series of sensors for measuring and acquiring an initial set of empirical 
values at various stages of the technical process while the technical process is 
5 operating based on a predetermined set of parameters; and 

a central processing unit being supplied with the initial set of empirical 
values, which screens a set of empirical values out of the initial set of empirical 
values using a screening algorithm to obtain a screened set of empirical values by: 
a division of the first set of empirical values in to classes based on a predetermined 

1 0 criteria, followed by an assessment of an empirical value within a class with respect 
to a predefined first threshold value, if a result of the assessment lies below a 
predefined first threshold value, then, the empirical value is screened out; and a 
further assessment of the class with respect to a predefined second threshold value; 
if a result of the further assessment lies below the second predefined threshold, then, 

15 screening out the class, 

said screened set of empirical values utilized for a simulation of the technical process 
to obtain a model result. 

REMARKS 

The foregoing amendments to the specification and claims under Article 41 
2 0 of the Patent Cooperation Treaty place the application into a form for prosecution 
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before the U.S. Patent and Trademark Office under 35 U.S.C §37L Accordingly, 
entry of these amendments before examination on the merits is hereby requested. 

Respectfully submitted, 



5 Steven H. Noll (reg. no. 28,982) 

Schiff Hardin & Waite 
Patent Department 
6600 Sears Tower 
Chicago, Illinois 60606 
10 Telephone: 312-258-5785 

CUSTOMER NO. 26574 



ATTORNEY FOR APPLICANT 
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SPECIFICATION 
TITLE 

METHOD AND ARRANGEMENT FOR REDUCING A NUMBER 
OF MEASURED VALUES OF A TECHNICAL SYSTEM 

5 BACKGROUND OF THE INVENTION 

Field of the Invention 

The present inventive method relates to a method and an arrangement for 
reducing a number of measured values of a technical system. 

m Description of the Related Art 

S10 A technical system, for example a process engineering plant, supplies a large 

w quantity of measured values per unit time, using different measured-value pickups 
;L (sensors). In the course of a number of days or weeks, a quantity of data accumulates 
Jlj which requires a correspondingly high computational power to process it. If the 
g measured values are to be used to adapt or reset the technical system, operations are 
15 often necessary whose complexity admits only some of the measured values. It is then 
a great disadvantage to select a specific number of measured values arbitrarily from the 
quantity of all measured values and to process them further, since measured values 
with low significance, for example measured values with a high degree of interference, 
have a considerable influence on the overall result and falsify the latter considerably. 
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SUMMARY OF THE INVENTION 

The object of the invention is to specify a method and an arrangement for 
reducing measured values, it being ensured that the remaining measured values have 
a high significance with regard to their description of the technical system. In order to 
achieve the object, a method of reducing a number of measured values of a technical 
system is specified, in which the measured values are divided into classes in 
accordance with predefined criteria. The measured values in a class are assessed and 
those measured values whose assessment lies below a predefined first threshold value 
are screened out. 

Screening out the measured values results in a reduction with regard to the 
number of measured values. A considerably reduced number of measured values is 
thus available for further processing. The further processing can be carried out with a 
lower computational outlay as compared with the non-reduced number of measured 
values. 

In one embodiment, the classes themselves are also assessed. In particular, a 
class whose assessment lies below a predefined second threshold value can be 
screened out. This results in an additional reduction in the number of measured values. 

In another embodiment, one criterion for the division into classes consists in that, 
for each class, measured values relating to a predefinition of setting parameters of the 
technical system are determined. Typically, the technical system is set using a 
predefined number of setting parameters and after the setting has been carried out, 
there follows a (mostly time-delayed) reaction of the system to the setting parameters 



(transient response, transient process of the system). After the setting operation, a 
specific quantity of measured values which can be associated with the transient process 
are therefore picked up, measured values which are associated with the predefined set 
of setting parameters continuing to accumulate after the transient process has 
concluded (transition to steady-state operation). By adjusting the setting parameters, 
a new class is determined. All the measured values which respectively accumulate 
after an adjustment to the setting parameters belong in their own class. Another 
development consists in that measured values within a class which can be associated 
with the respective transient process are screened out. In addition, erroneous 
measured values can be screened out. In many cases, the setting of large technical 
systems depends on long-term, steady-state operation. Measured values which relate 
to the transient process (of short duration in comparison with the steady-state operation 
after the transient process has concluded) are expediently screened out, since 
measured values for the steady-state operation are falsified by them. In particular 
within the context of modeling the technical system, the measured data about the 
steady-state behavior of the technical system are of interest. 

One refinement of this embodiment consists in reducing the number of 
measured values in a class in that at least one representative value for the measured 
values in the class is determined. Such a representative value may be: 

a) an average (e.g. a sliding average) of the measured values of the class, 

b) a maximum value of the measured values of the class, 

c) a minimum value of the measured values of the class, 



d) a median. 

In the case of variant (d), one advantage resides in the fact that it is always 
possible to determine a value which is actually present, while the average (a) does not 
itself occur as a value. 

Depending on the application, a suitable selection for determining the 
representative value of a class may be made. 

An entire class with measured values can be screened out if said class contains 
less than a predefined number of measured values. 

Another refinement consists in that those measured values are screened out 
which differ from a predefinable value by more than a predefined threshold value. 

The predefinable value can be an average of all the measured values of the 
class or a measured value to be expected in response to the respective setting 
parameters of the technical system. 

A development consists in that the measured values which are reduced in 
number are used for the simulation and/or the draft design of this or another technical 
system. 

The draft design of the technical system can in this case comprise both 
adaptation and redrafting of this or another technical system. 

In addition, within the context of the simulation, the behavior of the technical 
system can be checked by using the reduced number of measured values, with the 
objective of altering the system for developing a new system with altered predefinitions. 
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In addition, in order to achieve the object, an arrangement for reducing a number 
of measured values of a technical system is specified which has a processor unit, which 
processor unit is set up in such a way that the measured values can be divided into 
classes in accordance with predefined criteria. Measured values in a class can be 
5 assessed, measured values whose assessment lies below a predefined first threshold 
value are screened out. 

This arrangement is particularly suitable for carrying out the method according 
to the invention or one of its developments explained above. 

fi BRIEF DESCRIPTION OF THE DRAWINGS 

JO Fig- 1 shows a block diagram which contains steps in a method for reducing 

M : measured values. 

O Fig. 2 shows a schematic sketch of a recovery boiler. 

S Figs. 3 through 5 show input variables, actuating variables and outward variables 

of the recovery boiler. 



15 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 shows a block diagram which contains steps in a method of reducing 
measured values. In a first step 101 , the measured values are divided into classes. 
In particular, the division into classes is carried out with the effect that all the measured 
values which belong to one set of setting parameters are combined in one class. To 
20 this extent, each alteration to the setting parameters of the technical system found a 
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new class. In particular, altering the setting parameters is associated with a transient 
process of the technical system, this process, as opposed to a steady-state behavior, 
containing extreme fluctuations in the measured values. In a step 102, individual 
measured values within one class are screened out. These may be, for example, 
erroneous measured values, that is to say measured values which exhibit a high 
deviation with respect to the other measured values or an average of the measured 
values or measured values from the transient process. There are a number of possible 
ways of screening out individual measured values within a class: 

1 . Measured values which are too poor (based on a predefined comparative 
value); 

2. Measured values which belong to a transient process; 

3. Measured values which are erroneous; 

4. Determining a representative measured value as a representative for a 
plurality of measured values, in that the representative measured value 
is determined as an average of the measured values in a class or as a 
maximum value or a minimum value of these measured values. 

Measured values of this type are preferably not taken into account; they are 
removed from the respective class. This results in a considerable reduction in the 
number of measured values. In a step 1 03, individual classes are screened out. One 
criterion for screening out an entire class consists in that the class contains less than 
a predefined number of measured values. In a step 104, the measured values reduced 



in number are used for further processing. Further processing is, in particular, a 
simulation and/or a draft design of the technical system. 

Fig. 2 shows a schematic sketch of a recovery boiler. In the following text, by 
using the example of a "recovery boiler", an exemplary embodiment of the method 
described above will be illustrated. 

In the paper and pulp industry, various chemicals and also heat and electrical 
power are needed for the digestion of pulp. With the aid of the recovery boiler, the 
chemicals used and additional thermal energy may be recovered from a thickened 
waste process liquor (black liquor). A measure of the recovery of the chemicals is of 
critical importance for the economy of the overall plant. 

The black liquor is burned in a char bed 201 . In the process, an alkaline melt is 
produced and flows away via a line 202. In further process steps, the chemicals used 
are recovered from the constituents of the alkaline melt. Heat of combustion which is 
released is used to generate steam. The combustion of the waste liquor and therefore 
the recovery of the chemicals begins with the atomization of the black liquor via 
atomizer nozzles 204 into a combustion chamber 203. As they fall through the hot flue 
gas, particles of the atomized black liquor are dried. The dried liquor particles fall onto 
the char bed 201, first combustion and chemical reduction taking place. Volatile 
constituents and reaction products pass into an oxidation zone, in which oxidizing 
reactions proceed and in which the combustion is completed. 



Important objectives for the control of the recovery boiler are the steam 
production in order to obtain power, compliance with emission values from 
environmental points of view and the efficiency of the chemical reduction. 

The combustion operation, and therefore the objectives, are controlled in 
particular by the supply of air at three levels (primary air (PA), secondary air (SA), 
tertiary air (TA)). The overall process is subject to numerous influences, which have to 
be taken into account during the modeling: 

a) the measurement of the variables is subject to fluctuations which are 
often extreme; 

b) influencing variables which are not measured and cannot be measured 
exist; 

c) at each alteration to the settable parameters transient processes occur; 

d) the technical plant becomes soiled and is cleaned at predefined intervals, 
which has the effect of a drift over time in each case in the system 
behavior. 

The measured variables of the overall process are subdivided into input variables 
(cf. Fig. 3) and output variables (cf. Fig. 5). Measured variables are stored every 
minute. Four of the input variables are simultaneously also actuating variables (also: 
settable parameters; cf. Fig. 4). The actuating variables are to be viewed substantially 
as free parameters of the overall process which can be set independently of one 
another. Some of the other input variables are more or less dependent on the actuating 
variables. According to one predefinition, the variables "BL Front Pressure" and "BL 
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Back Pressure" are always to be regulated equally in the recovery boiler. The four 
actuating variables (cf. Fig. 4) are preferably to be stored as actuating variables (with 
the desired, preset value) and as input variables (with the measured, actual value). 

In the recovery boiler, one problem consists in the fact that, depending on the 
settable parameters, specific objectives, which are defined via measured variables, 
have to be met. Here, a three-stage procedure is selected in order to solve the 
problem: 

1 . The objectives to be considered are modeled by means of stochastic methods, 
these models being updated by means of new measurements (data-driven, 
empirical modeling). In this case, it is expedient to use not just a single model 
but global models for the identification of interesting areas in a parameter space 
determined by the objectives, and to use local models for the exact calculation 
of optimum operating points. The models used are assessed by means of 
quality measures. 

2. If the models considered are not sufficiently accurate because of the state of the 
data (quality measure), new operating points are deliberately evaluated in order 
to improve the model (experimental design). In addition, by using global 
stochastic optimization methods with regard to the objectives, attractive regions 
are identified on the basis of the current global model. 



3. For the local optimization, local models are constructed, and the data sets which 
are available are, if appropriate, deliberately expanded (experimental design). 
The objectives are physical/technical or economic criteria which, as a rule, have 
to meet boundary conditions and/or safety conditions. It is often the case that a number 
5 of these criteria have to be considered at the same time. A stochastic model can be 
used in particular for the purpose of simulating the objectives to be optimized and their 
dependence on the parameters to be set in the computer. This is necessary when 
measurements are very costly or very time-consuming. In the case of safety 
*S requirements, possible hazardous situations can be avoided. 

tjjo In the case of the recovery boiler, on-line optimization, which is based on a 

=p plurality of items of data, is necessary, since the physical/chemical processes cannot 
M be modeled quantitatively with sufficient accuracy and because the behavior of the 
y plant is subject to fluctuations in the course of operation. The knowledge about this 
behavior must continually be expanded by means of the deliberate selection of new 
r|5 operating points. Therefore, within the context of on-line optimization, the above- 
described three-stage procedure of stochastic modeling at mathematical optimization 
is to be recommended. 

DESCRIPTION OF THE INPUT VARIABLES 

The a input variables (a e N, N: set of natural numbers) generally depend on n 
20 actuating variables n e N and on random effects. They can be described as follows: 
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Let (Q, S, P) be a probability space and B v be a Borel a-algebra over R v (R: set of real 
numbers) for each v e N. The input variables are represented by a projection cp which 
can be measured via B n x S - B a : 

cp:R n xQ-R a (1 ). The definition set of the projection 

5 cp is a cartesian product of two sets. If one considers the respective projections onto 
the individual sets, then the following projections are obtained: 

(p x :Q — > R a , co — > <p(x, co) for all x e R n (2), 
y®:R n R a , x -> cp(x, co) for all co e Q (3) . 



= { } is a stochastic process having an index set R n and a projection cp w is 

]S a path in this stochastic process for each event wed. 

Qfc In the case of the recovery boiler, n=4 and a=14 (following the elimination of the 

variable "BL Back Pressure"). 

Because of the required ability to measure the projection cp x for each x e R n , the 

projection cp x is a random variable. Under suitable additional preconditions, expected 

values and higher moments can be considered. This access makes the step possible 
15 from stochastic models to deterministic optimization problems. In the case of a 

deterministic optimization problem, the target function can be set directly by means of 
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a variable, while the stochastic variable influences the target function but does not 
permit any deliberate setting. 

DESCRIPTION OF THE OUTPUT VARIABLES 

The process model M of the recovery boiler will be described as a function 
5 depending on the input variables and further random effects. In this case, let (Q, S, P) 
be the above probability space. The process model M is then a projection which can 
be 

;~ measured by B a x S - B b : 

f\ M : R a x fl -+ R b < 4 >< 

HO where b refers to the number of output variables. 

C Since the recovery boiler is subject to a cyclic drift over time (from cleaning 

s =i; phase to cleaning phase), a description using a time parameter is also conceivable. 
[2 The output variables may be represented by projections i that can be measured by B n 
x S - B b : 

V : R n x Q -> R b (5) , 

(x, ©) -> M(<p(x, cd), co) (6) . 
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If the respective projections onto the individual sets of the definition set are 
considered, then the following projections are obtained 

\j/ x : Q R b , o -> i|/(x, co) for all x e R n (7), 
M^iR 11 -> R b , x -> \|/( X/ co) for all CO € O (8) . 



3 v|/ x ; x e R n J 

q } is a stochastic process having an index set R n , and the projection qj w 

i=5 is a path in this stochastic process for each weQ. 
^ In the recovery boiler, b=15. 

3 The fact that, when defining ip, no distinction is drawn between the events oo 

used, does not mean that there is any restriction, since CI can be represented as a 
=7 cartesian product of an Q1 and an Q2. The above representation therefore also 
10 comprises the model: 

v)/ : R n x ft x x n 2 R b (9) ' 

(x, 01,02) -» 01), 02) (10). 
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DESCRIPTION OF THE DATA SETS AVAILABLE 

Using the descriptions in the two preceding sections, it is possible to combine the 
input variables and the output variables together to form measured variables 0. <t> is 
a projection that can be measured by B n x S - B m , where m = a + b, and 



c& : R n x CI -> R 



m 



(11) , 



(x, G>) 



V\|/(x, cd)J 



:i2) 



If the respective projections onto the individual sets of the definition set are 
considered again, then the following projections are obtained: 

<t> x :Q -» R m , co <E>(x, co) for all x e R n (13), 
^:R n -> R m , x -> <D(x, oo) for all co E fl (14) . 

(d> x ; x e R n ] 

{ } is a stochastic process with an index set R n and the projection O w is a 

path in this stochastic process for each co e Q. 
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For each chosen tuple x of actuating variables, a large number of 
implementations of O x in the recovery boiler are determined and stored, that is to say 
for each x } e R n , numerous implementations 

cp jk := 4>(xj,co jk ) {15) 
with ca j k € ; k = 1,2, . . * , v j ; 

v j € N ; j = 1,2, . . . , u ; u € H 



are considered. The stored data sets Djk of the recovery boiler are therefore (n + m) 
tuples: 

Djk = [o^J ' k = X ' 2 ' • • • ' v j 3 = l ' 2 > • ■ ■ ' u (16) • 

Here, D j1k1 is stored before D j2k2 if 
fal < 32) v fcl = 32) A ( k l < k 2)) 

DATA COMPRESSION BY DIVIDING THE CLASSES OF PARAMETERS 

Since, for each tuple x of actuating variables, there are generally a number of 
implementations of <t> x , because of the complex stochastic properties of the process to 
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be considered, the first step in the statistical data analysis is obviously to divide the 
classes of parameters by forming arithmetic averages. In addition, obviously erroneous 
data sets are separated out. An obviously erroneous data set is, for example, a 
physically impossible measurement which cannot possibly occur in real terms, in 
particular on the basis of a setting which has been made. 
Procedure: 

1 . Data sets in which the variable "BL Front Pressure" is not equal to the variable 
"BL Back Pressure" are screened out, since these two values must be equal 
according to the predefinition of the plant control system. The loss of data is 
very small. 

2. The data sets are divided into classes in which the four setting parameters (PA, 
SA, TA, BL Front Pressure, see above) are successively constant overtime, that 
is to say the jth class consists of the data sets Dj*. 

3. Classes in which there are fewer than 30 data sets are screened out, in order 
that transient processes do not have any great influence. 

4. For each class, an arithmetic average 0 j and an empirical standard deviation 
sj are determined for all the measured variables: 



(17) , 



'3 k = i 



v J k = l 



S-i - 



:i8) 



k = l 



5. Classes in which the averages forthe variables PA, SA, TA or BL Front Pressure 
are too far removed from the corresponding setting parameters are screened 
out. In these classes, therefore, the setting values could not be reached. 

CHARACTERISTIC STATISTICAL VARIABLES FOR THE GIVEN CLASSES AND 
THEIR GRAPHICAL REPRESENTATION 
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In addition to the arithmetic averages and the empirical standard deviations 
which have been determined for the individual classes, a common standard deviation 
s is further determined in accordance with 



s = 



j = i 



u 



3 = 1 



^2 



(19) 



£ 5 Here, u stands for the number of classes (205 here) and v for the sum of Vj, that 

^ is to say v is the number of all the measured values used (38,91 5 here). 

Although other modifications and changes may be suggested by those skilled 
in the art, it is the intention of the inventors to embody within the patent warranted 
hereon all changes and modifications as reasonably and properly come within the 
10 scope of their contribution to the art. 
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Abstract of the Disclosure 
A method and system for screening out empirically data collected from various 
steps of a technical process operating based on a first set of parameters, by utilizing 
a screening algorithm to reduce the size of the empirical data set, thus improving the 
modeling and revising of the technical process. The algorithm utilizes various classes 
associated with the empirical values and within each class performs an assessment 
with respect to preselected threshold values. The algorithm also performs an 
assessment with respect to another preselected threshold value, for the class as a 
whole. 
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Description SPECIFICATION 

Method and arrange ment for reducing a number of measured values of a 
t echnical s ystem TITLE 

T-he -METHOD AND ARRANGEMENT FOR REDUCING A NUMBER 
5 OF MEASURED VALUES OF A TECHNICAL SYSTEM 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present inventive method relates to a method and an arrangement for 

reducing a number of measured values of a technical system. 

yfO Description of the Related Art 

£ A technical system, for example a process engineering plant, supplies a large 

M- quantity of measured values per unit time, using different measured-value pickups 
JrJ (sensors). In the course of a number of days or weeks, a quantity of data accumulates 
J which requires a correspondingly high computational power to process it. If the 
j|5 measured values are to be used to adapt or reset the technical system, operations are 
often necessary whose complexity admits only some of the measured values. It is then 
a great disadvantage to select a specific number of measured values arbitrarily from the 
quantity of all measured values and to process them further, since measured values 
with low significance, for example measured values with a high degree of interference, 
20 have a considerable influence on the overall result and falsify the latter considerably. 
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SUMMARY OF THE INVENTION 

The object of the invention is to specify a method and an arrangement for 
reducing measured values, it being ensured that the remaining measured values have 
a high significance with regard to their description of the technical system. 

This object is achieved in accordance with the features of the independent patent 
claims. Deve l opments of the invention emerge from the dependent claims. In order to 
achieve the object, a method of reducing a number of measured values of a technical 
system is specified, inwhich in which the measured values are divided into classes in 
accordance with predefined criteria. The measured values in a class are assessed and 
those measured values whose assessment lies below a predefined first threshold value 
are screened out. 

Screening out the measured values results in a reduction with regard to the 
number of measured values. A considerably reduced number of measured values is 
thus available for further processing. The further processing can be carried out with a 
lower computational outlay as compared with the non-reduced number of measured 
values. 

One development consists in that In one embodiment the classes themselves 
are also assessed. In particular, a class whose assessment lies below a predefined 
second threshold value can be screened out. This results in an additional reduction in 
the number of measured values. 



Another development cons i sts in that In another embodiment one criterion for 
the division into classes consists in that, for each class, measured values relating to a 
predefinition of setting parameters of the technical system are determined. Typically, 
the technical system is set using a predefined number of setting parameters and after 
5 the setting has been carried out, there follows a (mostly time-delayed) reaction of the 
system to the setting parameters (transient response, transient process of the system). 
After the setting operation, a specific quantity of measured values which can be 
associated with the transient process are therefore picked up, measured values which 
y are associated with the predefined set of setting parameters continuing to accumulate 
flO after the transient process has concluded (transition to steady-state operation). By 
adjusting the setting parameters, a new class is determined. All the measured values 
M: which respectively accumulate after an 

O adjustment to the setting parameters belong in their own class. Another class. Another 
^ development consists in that measured values within a class which can be associated 
35 with the respective transient process are screened out. In addition, erroneous 
measured values can be screened out. In many cases, the setting of large technical 
systems depends on long-term, steady-state operation. Measured values which relate 
to the transient process (of short duration in comparison with the steady-state operation 
after the transient process has concluded) are expediently screened out, since 
20 measured values for the steady-state operation are falsified by them. In particular 
within the context of modeling the technical system, the measured data about the 
steady-state behavior of the technical system are of interest. 



One refinement of this embodiment consists in reducing the number of 
measured values in a class in that at least one representative value for the measured 
values in the class is determined. Such a representative value may be: 

a) an average (e.g. a sliding average) of the measured values of the class, 

b) a maximum value of the measured values of the class, 

c) a minimum value of the measured values of the class, 

d) a median. 

In the case of variant £d), one advantage resides in the fact that it is always 
possible to determine a value which is actually present, while the average (a) does not 
itself occur as a value. 

Depending on the application, a suitable selection for determining the 
representative value of a class may be made. 

An entire class with measured values can be screened out if said class contains 
less than a predefined number of measured values. 

Another refinement consists in that those measured values are screened out 
which differ from a predefinable value by more than a predefined threshold value. The 
value. 

The predefinable value can be an average of all the measured values of the 

class or a measured value to be expected in response to the respective setting 
parameters of the technical system. 



A development consists in that the measured values which are reduced in 
number are used for the simulation and/or the draft design of this or another technical 
system. 

The draft design of the technical system can in this case comprise both 
5 adaptation and redrafting of this or another technical system. 

In addition, within the context of the simulation, the behavior of the technical 
system can be checked by using the reduced number of measured values, with the 
objective of altering the system for developing a new system with altered predefinitions. 
^ In addition, in orderto achieve the object, an arrangement for reducing a number 

i^O of measured values of a technical system is specified which has a processor unit, which 
_rj processor unit is set up in such a way that the measured values can be divided into 
H classes in accordance with predefined criteria. Measured values in a class can be 
□ assessed, measured values whose assessment lies below a predefined first threshold 
Q value are screened out. 

35 This arrangement is particularly suitable for carrying out the method according 

to the invention or one of its developments explained above. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Eig Exemplary embodiments of the i nvention will be illustrated and explained 

below using the drawing, in whichl~ig . 1 shows a block diagram which contains steps 
20 in a method for reducing measured values^ 

Fig. 2 shows a schematic sketch of a recovery boiler^ 
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Figs. 3 -through 5 show input variables, actuating variables and outward 
variables of the recovery boiler. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 shows a block diagram which contains steps in a method of reducing 
5 measured values. In a first step 101 , the measured values are divided into classes. 
In particular, the division into classes is carried out with the effect that all the measured 
values which belong to one set of setting parameters are combined in one class. To 
3 this extent, each alteration to the setting parameters of the technical system found a 
% new class. In particular, altering the setting parameters is associated with a transient 
Jo process of the technical system, this process, as opposed to a steady-state behavior, 
II containing extreme fluctuations in the measured values. In a step 102, individual 

□ measured values within one class are screened out. These may be, for example, 

□ erroneous measured values, that is to say measured values which exhibit a high 
=f. deviation with respect to the other measured values or an average of the measured 
15 values or measured values from the transient process. There are a number of possible 

ways of screening out individual measured values within a class: 

1 . Measured values which are too poor (based on a predefined comparative 
value); 

2. Measured values which belong to a transient process; 
20 3. Measured values which are erroneous; 
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4. Determining a representative measured value as a representative for a 
plurality of measured values, in that the representative measured value 
is determined as an average of the 

measured values in a class or as a maximum value or a minimum value of these 

measured values. Measured values. 

Measured values of this type are preferably not taken into account; they are 

removed from the respective class. This results in a considerable reduction in the 
number of measured values. In a step 103, individual classes are screened out. One 
criterion for screening out an entire class consists in that the class contains less than 
a predefined number of measured values. In a step 104, the measured values reduced 
in number are used for further processing. Further processing is, in particular, a 
simulation and/or a draft design of the technical system. 

Fig. 2 shows a schematic sketch of a recovery boiler. In the following text, by 
using the example of a "recovery boiler", an exemplary embodiment of the method 
described above will be illustrated. 

In the paper and pulp industry, various chemicals and also heat and electrical 
power are needed for the digestion of pulp. With the aid of the recovery boiler, the 
chemicals used and additional thermal energy may be recovered from a thickened 
waste process liquor (black liquor). A measure of the recovery of the chemicals is of 
critical importance for the economy of the overall plant. 

The black liquor is burned in a char bed 201 . In the process, an alkaline melt is 
produced and flows away via a line 202. In further process steps, the chemicals used 
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are recovered from the constituents of the alkaline melt. Heat of combustion which is 
released is used to generate steam. The combustion of the waste liquor and therefore 
the recovery of the chemicals begins with the atomization of the black liquor via 
atomizer nozzles 204 into a combustion chamber 203. As they fall through the hot flue 
5 gas, particles of the atomized black liquor are dried. The dried liquor particles fall onto 
the char bed 201 , first combustion andchemical and chemical reduction taking place. 
Volatile constituents and reaction products pass into an oxidation zone, in which 
oxidizing reactions proceed and in which the combustion is completed. 
Q Important objectives for the control of the recovery boiler are the steam 

fip production in order to obtain power, compliance with emission values from 
S environmental points of view and the efficiency of the chemical reduction. 
£1 The combustion operation, and therefore the objectives, are controlled in 

p particular by the supply of air at three levels (primary air (PA), secondary air (SA), 
P tertiary air (TA)). The overall process is subject to numerous influences, which have to 
Ob be taken into account during the modeling: 

a) the measurement of the variables is subject to fluctuations which are 
often extreme; 

b) influencing variables which are not measured and cannot be measured 
exist; 

20 c) at each alteration to the settable parameters transient processes occur; 
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d) the technical plant becomes soiled and is cleaned at predefined intervals, 
which has the effect of a drift over time in each case in the system 
behavior. 

The measured variables of the overall process are subdivided into input variables 
(cf. Fig. 3) and output variables (cf. Fig. 5). Measured variables are stored every 
minute. Four of the input variables are simultaneously also actuating variables (also: 
settable parameters; cf. Fig. 4). The actuating variables are to be viewed substantially 
as free parameters of the overall process which can be set independently of one 
another. Some of the other i nput other input variables are more or less dependent on 
the actuating variables. According to one predefinition, the variables "BL Front 
Pressure" and "BL Back Pressure" are always to be regulated equally in the recovery 
boiler. The four actuating variables (cf. Fig. 4) are preferably to be stored as actuating 
variables (with the desired, preset value) and as input variables (with the measured, 
actual value). 

In the recovery boiler, one problem consists in the fact that, depending on the 
settable parameters, specific objectives, which are defined via measured variables, 
have to be met Here, a three-stage procedure is selected in order to solve the 
problem: 

1 . The objectives to be considered are modeled by means of stochastic methods, 
these models being updated by means of new measurements (data-driven, 
empirical modeling). In this case, it is expedient to use not just a single model 
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but global models for the identification of interesting areas in a parameter space 
determined by the objectives, and to use local models for the exact calculation 
of optimum operating points. The models used are assessed by means of 
quality measures. 

5 2. If the models considered are not sufficiently accurate because of the state of the 
data (quality measure), new operating points are deliberately evaluated in order 
to improve the model (experimental design). In addition, by using global 
S stochastic optimization methods with regard to the objectives, attractive regions 

are identified on the basis of the current global model. 

llo 3. Forthe local optimization, local models are constructed, and the data sets which 
O are available 

Q are, if appropriate, deliberately expanded (experimental design).The design). 

J^f The objectives are physical/technical or economic criteria which, as a rule, have 

to meet boundary conditions and/or safety conditions. It is often the case that a number 
15 of these criteria have to be considered at the same time. A stochastic model can be 
used in particular forthe purpose of simulating the objectives to be optimized and their 
dependence on the parameters to be set in the computer. This is necessary when 
measurements are very costly or very time-consuming. In the case of safety 
requirements, possible hazardous situations can be avoided. 
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In the case of the recovery boiler, on-line optimization, which is based on a 
plurality of items of data, is necessary, since the physical/chemical processes cannot 
be modeled quantitatively with sufficient accuracy and because the behavior of the 
plant is subject to fluctuations in the course of operation. The knowledge about this 
behavior must continually be expanded by means of the deliberate selection of new 
operating points. Therefore, within the context of on-line optimization, the 
above-described three-stage procedure of stochastic modeling at mathematical 
optimization is to be recommended. 

DESCRIPTION OF THE INPUT VARIABLES 

The a input variables (a e N, N: set of natural numbers) generally depend on n 
actuating variables n e N and on random effects. They can be described as follows: 

Let (Q, S, P) be a probability space and B v be a Borel a-algebra over RV (R: set of 
real numbers) for each v e N. The input variables are represented by a projection cp 
which can be measured via B n x S - B a : 

cp: R n x Q - R a WrTh e(1), The definition set of the 

projection cp is a cartesian product of two sets. If one considers the respective 
projections onto the individual sets, then the following projections are obtained: 

<p x : Q R a , co -+ <p(x, co) for all x <= R n (2), 
<ptt :fi n ^ R a^ x _^ <p( X/G) ) for all & 6 Q (3) . 
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(<p x ; x € R n ) 

1 } is a stochastic process having an index set 

R n and a projection cp w is a path in this stochastic process for each event weQ. 

In the case of the recovery boiler, n=4 and a=14 (following the elimination of the 

variable "BL Back Pressure"). 
5 Because of the required ability to measure the projection cp x for each x e R n , the 

projection cp x is a random variable. Under suitable additional preconditions, expected 

values and higher moments can be considered. This access makes the step possible 
yg from stochastic models to deterministic optimization problems. In the case of a 
U1 deterministic optimization problem, the target function can be set directly by means of 
jf 0 a variable, while the stochastic variable influences the target function but does not 

permit any deliberate setting. 

J DESCRIPTION OF THE OUTPUT VARIABLES 

M> The process model M of the recovery boiler will be described as a function 

depending on the input variables and further random effects. In this case, let (Q, S, P) 
15 be the above probability space. The process model M is then a projection which can 
be 

measured by B a x S - B b := Bbt 

M : R a x n R b (4) , 
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where b refers to the number of output variables. 

Since the recovery boiler is subject to a cyclic drift over time (from cleaning 
phase to cleaning phase), a description using a time parameter is also conceivable. 
The output variables may be represented by projections i that can be measured by B n 
x S - B b : 

V : R n x H -> R b (5) , 

(x, co) M(<p(x, co), co) (6) . 

If the respective projections onto the individual sets of the definition set are 
considered, then the following projections are obtained 

y x : n -> R b , to -> \\j(x,o>) for all x e R n (7), 
V^iR 11 -> R b , x v]/(x, co) for all co e £1 (8) . 



set R n , and the projection tt^ip^ is a path in this stochastic process for each ooefi. 




is a stochastic process having an index 
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In the recovery boiler, b=15. 

The fact that, when defining hjj, no distinction is drawn between the events u) 
used, does not mean that there is any restriction, since Q can be represented as a 
cartesian product of an Q1 and an Q2. The above representation therefore also 
comprises the model: 

(x, ©i, ©2) ©1), ©2) (10) • 
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DESCRIPTION OF THE DATA SETS AVAILABLE 

Using the descriptions in the two preceding sections, it is possible to combine the 
input variables and the output variables together to form measured variables O. O is 
a projection that can be measured by B n x S - B m , where m = a + b, and 



O : R n x CI R 



m 



:id , 



(x, co) 



^cp(x, co) 



(12) . 



If the respective projections onto the individual sets of the definition set are 
considered again, then the following projections are obtained: 

3> x : Q ~> R m , co cp(x, co) for all x e R n (13), 
O^rR 11 -> R m , x 3>(x,co) for all co e Q (14). 

{a> x ; x g R n \ 

1 ; is a stochastic process with an index set R n 

and the projection O w is a path in this stochastic process for each co e Q. 
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For each chosen tuple x of actuating variables, a large number of 
implementations of <J> X in the recovery boiler are determined and stored, that is to say 
for each Xj e R n , numerous implementations 

(D jk := *(xj,©j k ) (15) 
with cd j € Q ; k = 1,2, . . . , Vj ; 

v j e N ; j = 1,2, — , u ; u e S 



are considered. The stored data sets Djk of the recovery boiler are therefore (n + m) 
tuples: 



D- 



jk 



k = 1,2, 



v- 



j = 1,2, . - . , u 



(16) 



Here, D j1k1 is stored before D j2k2 if 
(di < J2) v (61 = 12) A ( k i < ^2)) 

DATA COMPRESSION BY DIVIDING THE CLASSES OF PARAMETERS 

Since, for each tuple x of actuating variables, there are generally a number of 
implementations of O x , because of the complex stochastic properties of the process to 
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be considered, the first step in the statistical data analysis is obviously to divide the 
classes of parameters by forming arithmetic averages. In addition, obviously erroneous 
data sets are separated out. An obviously erroneous data set is, for example, a 
physically impossible measurement which cannot possibly occur in real terms, in 
particular on the basis of a setting which has been made. 
Procedure: 

1 . Data sets in which the variable "BL Front Pressure" is not equal to the variable 
"BL Back Pressure" are screened out, since these two values must be equal 
according to the predefinition of the plant control system. The loss of data is 
very small. 

2. The data sets are divided into classes in which the four setting parameters (PA, 
SA, TA, BL Front Pressure, see above) are successively constant overtime, that 
is to say the jth class consists of the data sets D *. 

3. Classes in which there are fewer than 30 data sets are screened out, in order 
that transient processes do not have any great influence. 

4. For each class, an arithmetic average O j and an empirical standard deviation 
sj are determined for all the measured variables: 
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5. Classes in which the averages for the variables PA, SA, TA or BL Front Pressure 
are too far removed from the corresponding setting parameters are screened 
out. In these classes, therefore, the setting values could not be 
5 rcachcd.CI IARACTCRIGTIC reached. 

CHARACTERISTIC STATISTICAL VARIABLES FOR THE GIVEN CLASSES AND 
THEIR GRAPHICAL REPRESENTATION 
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In addition to the arithmetic averages and the empirical standard deviations 
which have been determined for the individual classes, a common standard deviation 
s is further determined in accordance with 



Here, u stands for the number of classes (205 here) and v for the sum of Vj, that 
is to say v is the number of all the measured values used (38,91 5 here). 

Abstract Although other modifications and changes may be su gg es ted b y 
those skilled in the art, it is the intention of the invento rs to embody within the 
patent warranted hereon all changes and modifications as reasonably and 
properly come within the scope of their contribut ion to the art. 




(19) 



s = 



1 
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M e thod and arrangement for reducing a number of measured values of a technical 
system Abstract of the Disclosure 

A method of reducing a number of measured va l ues of the techn i cal system is 
specified in which the measured values are divided into classes i n accordance with 
predefined criteria. The measured va l ues in a class are assessed and those measured 
values whose assessment lies below a predefined f i rst threshold value are screened 
otrt A method and system for screening out empirically data collected from 
various steps of a technical process operating based on a first set of parameters, 
by utilizing a screening algorithm to reduce the size of the empirical data set 
thus improving the modeling and revising of the technical process. The 
al gorithm utilizes various classes associated with the empirical values and within 
each class performs an assessment with respect to preselected threshold values- 
The algorithm also performs an assessment with respect to another preselected 
threshold value, for the class as a whole . 
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Method and arrangement for reducing a number of 
measured values of a technical system 

The method relates to a method and an arrangement for 
reducing a number of measured values of a technical 
system. 



A technical system, for example a process engineering 
plant, supplies a large quantity of measured values per 
unit time, using different measured-value pickups 
(sensors) . In the course of a number of days or weeks, 
a quantity of data accumulates which requires a 
correspondingly high computational power to process it. 
If the measured values are to be used to adapt or reset 
the technical system, operations are often necessary 
whose complexity admits only some of the measured 
values. It is then a great disadvantage to select a 
specific number of measured values arbitrarily from the 
quantity of all measured values and to process them 
further, since measured values with low significance, 
for example measured values with a high degree of 
interference, have a considerable influence on the 
overall result and falsify the latter considerably. 

The object of the invention is to specify a method and 
an arrangement for reducing measured values, it being 
ensured that the remaining measured values have a high 
significance with regard to their description of the 
technical system. 

This object is achieved in accordance with the features 
of the independent patent claims. Developments of the 
invention emerge from the dependent claims. 
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in order to achieve the object, a method of reducing 
number of measured values of a technical system 
specified, in 
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which the measured values are divided into classes in 
accordance with predefined criteria. The measured 
values in a class are assessed and those measured 
values whose assessment lies below a predefined first 
threshold value are screened out. 

Screening out the measured values results in a 
reduction with regard to the number of measured values. 
A considerably reduced number of measured values is 
thus available for further processing. The further 
processing can be carried out with a lower 
computational outlay as compared with the non-reduced 
number of measured values. 



One development consists in that the classes themselves 
are also assessed. In particular, a class whose 
assessment lies below a predefined second threshold 
value can be screened out. This results in an 
additional reduction in the number of measured values. 

Another development consists in that one criterion for 
the division into classes consists in that, for each 
class f measured values relating to a predef inition of 
setting parameters of the technical system are 
determined. Typically, the technical system is set 
using a predefined number of setting parameters and 
after the setting has been carried out, there follows a 
(mostly time-delayed) reaction of the system to the 
setting parameters (transient response, transient 
process of the system) . After the setting operation, a 
specific quantity of measured values which can be 
associated with the transient process are therefore 
picked up, measured values which are associated with 
the predefined set of setting parameters continuing to 
accumulate after the transient process has concluded 
(transition to steady-state operation) . By adjusting 
the setting parameters, a new class is determined. All 
the measured values which respectively accumulate after 
an 



adjustment to the setting parameters belong in their 
own class. 
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Another development consists in that measured values 
within a class which can be associated with the 
respective transient process are screened out. In 
addition, erroneous measured values can be screened 
out. In many cases, the setting of large technical 
systems depends on long-term, steady-state operation. 
Measured values which relate to the transient process 
(of short duration in comparison with the steady-state 
operation after the transient process has concluded) 
are expediently screened out, since measured values for 
the steady-state operation are falsified by them. In 
particular within the context of modeling the technical 
system, the measured data about the steady-state 
behavior of the technical system are of interest. 

One refinement consists in reducing the number of 
measured values in a class in that at least one 
representative value for the measured values in the 
class is determined. Such a representative value may 
be: 

a) an average (e.g. a sliding average) of the 
measured values of the class, 

b) a maximum value of the measured values of the 
class , 

c) a minimum value of the measured values of the 
class , 

d) a median. 

In the case of variant d) , one advantage resides in the 
fact that it is always possible to determine a value 
which is actually present, while the average a) does 
not itself occur as a value. 

Depending on the application, a suitable selection for 
determining the representative value of a class may be 
made . 
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An entire class with measured values can be screened 
out if said class contains less than a predefined 
number of measured values. 

Another refinement consists in that those measured 
values are screened out which differ from a 
predefinable value by more than a predefined threshold 
value . 
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The predefinable value can be an average of all the 
measured values of the class or a measured value to be 
expected in response to the respective setting 
parameters of the technical system. 

A development consists in that the measured values 
which are reduced in number are used for the simulation 
and/or the draft design of this or another technical 
system. 

The draft design of the technical system can in this 
case comprise both adaptation and redrafting of this or 
another technical system. 

In addition, within the context of the simulation, the 
behavior of the technical system can be checked by 
using the reduced number of measured values, with the 
objective of altering the system for developing a new 
system with altered predef initions . 

In addition, in order to achieve the object, an 
arrangement for reducing a number of measured values of 
a technical system is specified which has a processor 
unit, which processor unit is set up in such a way that 
the measured values can be divided into classes in 
accordance with predefined criteria. Measured values in 
a class can be assessed, measured values whose 
assessment lies below a predefined first threshold 
value are screened out. 

This arrangement is particularly suitable for carrying 
out the method according to the invention or one of its 
developments explained above. 

Exemplary embodiments of the invention will be 
illustrated and explained below using the drawing, in 
which 
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Fig. 1 shows a block diagram which contains steps in 
a method for reducing measured values ; 

Fig. 2 shows a schematic sketch of a recovery 
boiler; 

Figs. 3-5 show input variables, actuating variables and 
outward variables of the recovery boiler. 

Fig. 1 shows a block diagram which contains steps in a 
method of reducing measured values. In a first step 
101, the measured values are divided into classes. In 
particular, the division into classes is carried out 
with the effect that all the measured values which 
belong to one set of setting parameters are combined in 
one class. To this extent, each alteration to the 
setting parameters of the technical system found a new 
class. In particular, altering the setting parameters 
is associated with a transient process of the technical 
system, this process, as opposed to a steady-state 
behavior, containing extreme fluctuations in the 
measured values. In a step 102, individual measured 
values within one class are screened out. These may be, 
for example, erroneous measured values, that is to say 
measured values which exhibit a high deviation with 
respect to the other measured values or an average of 
the measured values or measured values from the 
transient process. There are a number of possible ways 
of screening out individual measured values within a 
class : 

1. Measured values which are too poor (based on a 
predefined comparative value) ; 

2. Measured values which belong to a transient 
process; 

3. Measured values which are erroneous; 

4. Determining a representative measured value as 
a representative for a plurality of measured 
values, in that the representative measured 
value is determined as an average of the 
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measured values in a class or as a maximum 
value or a minimum value of these measured 
values . 



GR 98 P 5866 

- 6 - 

Measured values of this type are preferably not taken 
into account; they are removed from the respective 
class. This results in a considerable reduction in the 
number of measured values. In a step 103, individual 
classes are screened out. One criterion for screening 
out an entire class consists in that the class contains 
less than a predefined number of measured values. In a 
step 104, the measured values reduced in number are 
used for further processing. Further processing is, in 
particular, a simulation and/or a draft design of the 
technical system. 

Fig. 2 shows a schematic sketch of a recovery boiler. 
In the following text, by using the example of a 
"recovery boiler", an exemplary embodiment of the 
method described above will be illustrated. 

In the paper and pulp industry, various chemicals and 
also heat and electrical power are needed for the 
digestion of pulp. With the aid of the recovery boiler, 
the chemicals used and additional thermal energy may be 
recovered from a thickened waste process liquor (black 
liquor) . A measure of the recovery of the chemicals is 
of critical importance for the economy of the overall 
plant . 

The black liquor is burned in a char bed 201. In the 
process, an alkaline melt is produced and flows away 
via a line 202. In further process steps, the chemicals 
used are recovered from the constituents of the 
alkaline melt. Heat of combustion which is released is 
used to generate steam. The combustion of the waste 
liquor and therefore the recovery of the chemicals 
begins with the atomization of the black liquor via 
atomizer nozzles 204 into a combustion chamber 203. As 
they fall through the hot flue gas, particles of the 
atomized black liquor are dried. The dried liquor 
particles fall onto the char bed 201, first combustion 
and 
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chemical reduction taking place. Volatile constituents 
and reaction products pass into an oxidation zone, in 
which oxidizing reactions proceed and in which the 
combustion is completed. 

Important objectives for the control of the recovery 
boiler are the steam production in order to obtain 
power, compliance with emission values from 
environmental points of view and the efficiency of the 
chemical reduction . 

The combustion operation, and therefore the objectives, 
are controlled in particular by the supply of air at 
three levels (primary air (PA), secondary air (SA) , 
tertiary air (TA) ) . The overall process is subject to 
numerous influences, which have to be taken into 
account during the modeling: 

a) the measurement of the variables is subject to 
fluctuations which are often extreme; 

b) influencing variables which are not measured 
and cannot be measured exist ; 

c) at each alteration to the settable parameters 
transient processes occur; 

d) the technical plant becomes soiled and is 
cleaned at predefined intervals, which has the 
effect of a drift over time in each case in the 
system behavior. 

The measured variables of the overall process are 
subdivided into input variables (cf. Fig. 3) and output 
variables (cf . Fig. 5) . Measured variables are stored 
every minute. Four of the input variables are 
simultaneously also actuating variables (also: settable 
parameters; cf . Fig. 4) . The actuating variables are to 
be viewed substantially as free parameters of the 
overall process which can be set independently of one 
another. Some of the other 
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input variables are more or less dependent on the 
actuating variables. According to one predef inition, 
the variables "BL Front Pressure" and M BL Back 
Pressure" are always to be regulated equally in the 
recovery boiler. The four actuating variables (cf. Fig. 
4) are preferably to be stored as actuating variables 
(with the desired, preset value) and as input variables 
(with the measured, actual value) . 

In the recovery boiler, one problem consists in the 
fact that, depending on the settable parameters, 
specific objectives, which are defined via measured 
variables, have to be met. Here, a three-stage 
procedure is selected in order to solve the problem: 

1. The objectives to be considered are modeled by 
means of stochastic methods, these models being 
updated by means of new measurements (data-driven, 
empirical modeling) . In this case, it is expedient 
to use not just a single model but global models 
for the identification of interesting areas in a 
parameter space determined by the objectives, and 
to use local models for the exact calculation of 
optimum operating points. The models used are 
assessed by means of quality measures. 

2. If the models considered are not sufficiently 
accurate because of the state of the data (quality 
measure) , new operating points are deliberately 
evaluated in order to improve the model 
(experimental design) . In addition, by using 
global stochastic optimization methods with regard 
to the objectives, attractive regions are 
identified on the basis of the current global 
model . 

3. For the local optimization, local models are 
constructed, and the data sets which are available 
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are, if appropriate, 
(experimental design) . 



deliberately expanded 
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The objectives are physical/technical or economic 
criteria which, as a rule, have to meet boundary 
conditions and/or safety conditions. It is often the 
case that a number of these criteria have to be 
5 considered at the same time. A stochastic model can be 
used in particular for the purpose of simulating the 
objectives to be optimized and their dependence on the 
parameters to be set in the computer. This is necessary 
when measurements are very costly or very time- 
10 consuming. In the case of safety requirements, possible 
hazardous situations can be avoided. 

In the case of the recovery boiler, on-line 
optimization, which is based on a plurality of items of 

15 data, is necessary, since the physical/chemical 
processes cannot be modeled quantitatively with 
sufficient accuracy and because the behavior of the 
plant is subject to fluctuations in the course of 
operation. The knowledge about this behavior must 

20 continually be expanded by means of the deliberate 
selection of new operating points. Therefore, within 
the context of on-line optimization, the above- 
described three-stage procedure of stochastic modeling 
at mathematical optimization is to be recommended. 

25 

DESCRIPTION OF THE INPUT VARIABLES 

The a input variables (a e N, N: set of natural 
numbers) generally depend on n actuating variables n e 
30 N and on random effects. They can be described as 
follows : 

Let (Q, S, P) be a probability space and B v be a Borel 
o-algebra over R v (R: set of real numbers) for each v e 
35 N. The input variables are represented by a projection 
cp which can be measured via B 11 x S - B*: 



cp: R n x Q R a 



(1) 
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The definition set of the projection 9 is a cartesian 
product of two sets. If one considers the respective 
projections onto the individual sets, then the 
following projections are obtained: 

5 

cp x : Q -> R a , co -> (p(x, o)) for all x € R n (2), 
cp^iR 11 -» R a , x -> q>(x, ©) for all co efl (3) . 

tp x ; x e R n | is a stochastic process having an index set 

R n and a projection (p 01 is a path in this stochastic 
process for each event 0 e Q. 

10 

In the case of the recovery boiler, n=4 and a=14 
(following the elimination of the variable U BL Back 
Pressure") . 

15 Because of the required ability to measure the 
projection <p x for each x e R n , the projection <p x is a 
random variable. Under suitable additional 

preconditions, expected values and higher moments can 
be considered. This access makes the step possible from 

20 stochastic models to deterministic optimization 
problems. In the case of a deterministic optimization 
problem, the target function can be set directly by 
means of a variable, while the stochastic variable 
influences the target function but does not permit any 

25 deliberate setting. 

DESCRIPTION OF THE OUTPUT VARIABLES 

The process model M of the recovery boiler will be 
30 described as a function depending on the input 
variables and further random effects. In this case, let 
(Q, S, P) be the above probability space. The process 
model M is then a projection which can be measured by 
B* x S - B* 5 : 



35 
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M : R a x n R b (4) , 

where b refers to the number of output variables. 

Since the recovery boiler is subject to a cyclic drift 
over time (from cleaning phase to cleaning phase) , a 
description using a time parameter is also conceivable. 
The output variables may be represented by projections 
\y that can be measured by B n x S - B b : 



\\f : R n x fi -> R x 



(x, co) M(q>(x, cd), co) 



(5) , 

(6) . 



If the respective projections onto the individual sets 
of the definition set are considered, then the 
following projections are obtained 



V/ x : n -> R b , o> -> \\r(x, co) for all x e R n (7), 



i\^:R n -> R b , x -> m/(x, co) for all co € Q 



(8) 



v|/ x ; x e R | is a stochastic process having an index 

set R n , and the projection \\f° is a path in this 
stochastic process for each co e Q. 



In the recovery boiler, b=15. 



The fact that, when defining no distinction is drawn 
between the events co used, does not mean that there is 
any restriction, since Q can be represented as a 
cartesian product of an and an CI2 - The above 

representation therefore also comprises the model: 
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(x, 01,02) -» m(<p(x,Oi), o 2 ) 



(9), 
(10) 



DESCRIPTION OF THE DATA SETS AVAILABLE 



10 



Using the descriptions in the two preceding sections, 
it is possible to combine the input variables and the 
output variables together to form measured variables O. 
O is a projection that can be measured by B 11 x S - B" 1 , 
where m = a + b, and 



0) 



R n x Q 



(11) , 



{x, o) 



r cp(x, to) 
l\|/(x, ©), 



(12) 



If the respective projections onto the individual sets 
of the definition set are considered again, then the 
15 following projections are obtained: 

O x ; a -» R m , o ~> 0>(x, o) for all x e R n (13), 
O^rR 11 -> R m , x -» <D(x,o)) for all o e Q (14) . 

<D X ; x g R n | is a stochastic process with an index set R n 

and the projection O* 0 is a path in this stochastic 
20 process for each 00 e Q. 

For each chosen tuple x of actuating variables, a large 
number of implementations of O x in the recovery boiler 
are determined and stored, that is to say for each Xj e 
25 R n , numerous implementations 
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d> jk := 4>(xj, co jk ) 



(15) 



with coj^ € Q ; k = 1,2, 



v j e N ; j = 1,2, . . . , u ; u e B 



are considered. The stored data sets Dj k of the recovery 
5 boiler are therefore (n + m) tuples: 



DATA COMPRESSION BY DIVIDING THE CLASSES OF PARAMETERS 

15 Since, for each tuple x of actuating variables, there 
are generally a number of implementations of O x/ 
because of the complex stochastic properties of the 
process to be considered, the first step in the 
statistical data analysis is obviously to divide the 

20 classes of parameters by forming arithmetic averages. 
In addition, obviously erroneous data sets are 
separated out. An obviously erroneous data set is, for 
example, a physically impossible measurement which 
cannot possibly occur in real terms, in particular on 

25 the basis of a setting which has been made. 

Procedure : 





(16) . 



Here, Dj ik i is stored before Dj 2 k2 if 



10 



(ji < 32) v ((31 = j 2 ) a (ki < k 2 )) 



30 



1. Data sets in which the variable NX BL Front 
Pressure" is not equal to the variable "BL Back 
Pressure" are 
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screened out, since these two values must be equal 
according to the predef inition of the plant 
control system. The loss of data is very small. 



5 2, 



10 



The data sets are divided into classes in which 
the four setting parameters (PA, SA, TA, BL Front 
Pressure, see above) are successively constant 
over time, that is to say the jth class consists 
of the data sets Dj.. 

Classes in which there are fewer than 30 data sets 
are screened out, in order that transient 
processes do not have any great influence. 



15 4. For each class, an arithmetic average Oj and an 
empirical standard deviation sj are determined for 
all the measured variables: 



1 3 
3 k = l 



(17), 



3 k = l J 



k = l 



> J 



(18) . 



20 



25 



5. Classes in which the averages for the variables 
PA, SA, TA or BL Front Pressure are too far 
removed from the corresponding setting parameters 
are screened out. In these classes, therefore, the 
setting values could not be reached. 
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CHARACTERISTIC STATISTICAL VARIABLES FOR THE GIVEN 
CLASSES AND THEIR GRAPHICAL REPRESENTATION 

In addition to the arithmetic averages and the 
5 empirical standard deviations which have been 
determined for the individual classes, a common 
standard deviation s is further determined in 
accordance with 




s 



(19) 



1 




\ 



1 



v 



- 1 



J 



10 



Here, u stands for the number of classes (205 here) and 
v for the sum of Vj, that is to say v is the number of 
all the measured values used (38,915 here). 
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1. A method of reducing a number of measured 
variables of a technical system, 

a) in which the measured values are divided into 
classes in accordance with predefined criteria; 

b) in which the measured values in a class are 
assessed and measured values whose assessment 
lies below a predefined first threshold value 
are screened out; 

c) in which the classes are assessed and a class 
for which the assessment lies below a 
predefined second threshold value is screened 
out . 



2. The method as claimed in one of the preceding 
claims, in which one criterion for the division 
into classes consists in that, for each class, 
measured values relating to a predef inition of 
setting parameters of the technical system are 
determined. 



3. The method as claimed in one of the preceding 
claims, in which, in one class, measured values 
relating to a transient process and/or erroneous 
measured values are determined and screened out. 

4. The method as claimed in one of the preceding 
claims, in which, in a class, the number of 
measured values is reduced in that at least one 
representative value for the measured values in 
the class is determined. 

5. The method as claimed in claim 4, in which the 
representative value is determined as 

a) an average of the measured values in the class, 
or 

b) a maximum value or a minimum value of the 
measured values in the class; 
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6. The method as claimed in one of the preceding 
claims, in which a class which has fewer than a 
predefined number of measured values is screened 
out . 

5 

7 . The method as claimed in one of the preceding 
claims, in which, in a class, those measured 
values which differ from a predefinable value by 
more than a predefinable threshold value are 

10 screened out. 

8 . The method as claimed in one of the preceding 
claims, in which the reduced measured values are 
used for the simulation and/or for the draft 

15 design of the technical system. 

9. An arrangement for reducing a number of measured 
values of a technical system, 

having a processor unit which is set up in such a 
20 way that 

a) the measured values are divided into classes in 
accordance with predefined criteria; 

b) measured values in a class can be assessed and 
measured values whose assessment lies below a 

25 predefined first threshold value are screened 

out ; 

c) the classes are assessed and a class for which 
the assessment lies below a predefined second 
threshold value is screened out. 

30 

10. The arrangement as claimed in claim 9, in which 
the processor unit is set up in such a way that 
the classes are assessed and a class for which the 
assessment lies below a predefined second 

35 threshold value is screened out. 



AMENDED SHEET 
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Abstract 

Method and arrangement for reducing a number of 
measured values of a technical system 

A method of reducing a number of measured values of the 
technical system is specified in which the measured 
values are divided into classes in accordance with 
predefined criteria. The measured values in a class are 
assessed and those measured values whose assessment 
lies below a predefined first threshold value are 
screened out. 
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FIG1 



Dividing the measured values 
into classes 



Screening out individual measured 
values within a class 
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Screening out individual classes 



Further processing of the 
measured values 
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FIG 3 
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FIG 5 
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der Abwicklung aller damit verbundenen Geschafte 
vor dem Patent- und Warenzeichenamt: (Name und 
Registrationsnummer anfuhren) 



POWER OF ATTORNEY As a named inventor, I 
hereby appoint the following attorney(s) and/or 
agent(s) to prosecute this application and transact all 
business in the Patent and Trademark Office connec- 
ted therewith. (list name and registration number) 



And I hereby appoint 

Messrs. John D. Simpson (Registration No. 19.842 ) Lewis T. Steadman (17,074 ). William C Slueberj(16,453), P Phillips Connor X1&259), Dennts A Gross 
(-24-44©), Marvin Moody,£l&549), Steven H. Noil ( 28,982) , Brett A. Valiquet(2I£41), Thomas I. Ros s (29,2? 5), Kevin W. Guynn (29.927 ). Edward A Lehmann 
<22^42), James D. Hobart (24J49), Robert M. Barrett. (30, 142), James Van Santen (16,584 ), J. Arthur Gross (13,615), Richard J. Schwarz (13 47? ) and 
Melvin A. Robinson David R. Metzger.Q2,ft1-9), John R. Garrett_£2ZjaS8) all members of the firm of Hill, Steadman & Simpson, A Professional Corpo- 

ration. 



Telefongesprache bitte richten an: 
(Name und Teiefonnummer) 



Direct Telephone Calls to: (name and telephone num- 
ber) 

312/876-0200 

Ext 



Postanschrift: Send Correspondence to: 

HILL. STEADMAN & SIM PSON 
A Professional C orporation 
85th Floor Sears Tower, Chicago, Illinois L 606 06_ 



Voller Name des einzigen oder ursprunglichen Erfinders: 

SOHAFFI FR r Stefan 


Full name of sole or first inventor: 


U ntersc h rift des Erfinders Datum 

c/U- M M- M. }j 


Inventor's signature Date 


Wohns^z J 

D-861 99 A/jasbura. Germany JDS X. 


Residence 


Staatsangehbrigkeit 

Bundesrepublik Deutschland 


Citizenship 


Postanschrift 

Paul-Lincke-Str. 15 


Post Office Addess 


D-861 99 Augsburg 
Bundesrepublik Deutschland 




Voller Name des zweiten Miterfinders (falls zutreffend): 

STURM, Thomas 


Full name of second joint inventor, if any; 


Unterscht^kle^ Eifinde^ Datum 

<Mo MM 53 


Second Inventor's signature Date 


-ift/ohnsitz 

D-81673 Munchen, Germany T)(5 XL 


Residence 


Staatsangehorigkeit 

Bundesrepublik Deutschland 


Citizenship 


Postanschrift 

Schlusselbergstr. 16 


Post Office Address 


D-81673 Munchen 
Bundesrepublik Deutschland 





(Bitte entsprechende informationen und Unterschriften im 
Falfe von dritten und weiteren Miterfmdern angeben). 



(Supply similar information and signature for third and 
subsequent joint inventors). 
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